Optically mediated nonlinear quantum optomechanics 
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We consider theoretically the optomechanical interaction of several mechanical modes with a single 
quantized cavity field mode for linear and quadratic coupling. We focus specifically on situations 
where the optical dissipation is the dominant source of damping, in which case the optical field can 
be adiabatically eliminated, resulting in effective multimode interactions between the mechanical 
modes. In the case of linear coupling, the coherent contribution to the interaction can be exploited 
e.g. in quantum state swapping protocols, while the incoherent part leads to significant modifications 
of cold damping or amplification from the single-mode situation. Quadratic coupling can result 
in a wealth of possible effective interactions including the analogs of second-harmonic generation 
and four-wave mixing in nonlinear optics, with specific forms depending sensitively on the sign 
of the coupling. The cavity-mediated mechanical interaction of two modes is investigated in two 
limiting cases, the resolved sideband and the Doppler regime. As an illustrative application of the 
formal analysis we discuss in some detail a two-mode system where a Bose-Einstein condensate is 
optomechanically linearly coupled to the moving end mirror of a Fabry-Perot cavity. 



I. INTRODUCTION 

Quantum optomechanics has recently witnessed a 
rapid succession of key advances, with important mile- 
stones including the cooling of several systems to the 
quantum ground state of vibration the demon- 

stration of strong optomechanical coupling and of 

quantum coherent coupling between optical and phonon 
fields , the realization of optomechanical systems from 
ultracold atomic ensembles @, Q, and impressive ad- 
vances toward quantum state transfer [10l - ll5l |. to men- 
tion just a few examples. These developments hint at 
the promise for optomechanics applications in quantum 
metrology 0, [l7| , quantum information science fl8l - l2C| , 
and also in possible tests of fundamental physics ques- 
tions [ll|,[22|. So far, though, the focus of most advances 
has been on single mode dynamics of the oscillator [23l — 
l25j |. and only few authors have started to investigate 
the effects of multiple mechanical modes [2(| l26l - l28j . In 
complete analogy with the situation in quantum optics, 
however, single-mode quantum acoustics fails to capture 
many of its most important aspects, including most obvi- 
ously perhaps propagation phenomena, nonlinear acous- 
tical effects, focusing and defocusing of acoustic waves, 
etc. 

This paper represents a step toward the development of 
multimode quantum acoustics. The specific situation we 
consider is a multimode mechanical field - perhaps pro- 
vided by several mechanical oscillators or atomic samples 
trapped near extrema of an optical lattice, the mode cou- 
pling being provided by a common optical field. In some 
ways, this can be thought of as describing a situation 
analogous to collisions in atomic systems. Very much like 
collisions are effective interactions between atoms result- 
ing from the elimination of the electromagnetic field, the 
optically mediated effective coupling between phononic 
modes can also be thought of in these terms. Gener- 
ally speaking, collisions tend to be largely incoherent 
processes resulting e.g. in decoherence and dephasing in 
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FIG. 1. Non-absorptive dielectric membranes inside a Fabry- 
Perot resonator which interact with a single mode cavity field. 



atomic ensembles, but it is known from ultracold atomic 
systems that this needs not be the case. In ultracold 
systems, such as quantum degenerate atomic gases, col- 
lisions can be a coherent process, and have been demon- 
strated to lead to nonlinear atom optics effects [2£| |3(| 
such as e.g matter-wave four-wave mixing, soliton gener- 
ation, the creation of entangled atom pairs, and more. In 
analogy to that situation, one can expect similar effects 
to occur in quantum acoustics. 

Our main goal is to develop a simple theory that shows 
that this is indeed the case, and that incoherent effects 
due to quantum fluctuations - associated primarily with 
the elimination of the electromagnetic field - need not 
overwhelm coherent effects. This paves the way to in- 
vestigate nonlinear effects in phonon physics reminiscent 
of those encountered in nonlinear optics and nonlinear 
atom optics; with the promise of similarly exciting appli- 
cations. 

Our starting point is a relatively general discussion of 
optically mediated effective interactions between many 
mechanical modes, discussing in some detail the quan- 
tum noise and damping associated with that interaction, 
both for linear and quadratic optomechanical coupling. 
We identify a number of specific forms of the effective 



2 



interaction that lead to qualitatively different dynamics, 
in particular a regime where the coupling of the quan- 
tum fluctuations of the mechanical modes is intrinsically 
nonlinear. We then illustrate the general analysis with 
an example involving just two mechanical modes, with 
emphasis on two limiting cases, the resolved sideband 
and the Doppler regimes. A specific example is given by 
the state transfer between two mechanical modes, for ex- 
ample a Bose-Einstein condensate coupled to the center- 
of-mass mode of a moving end mirror of a Fabry-Perot 
resonator through a cavity field. We comment on the 
optimization of the transfer fidelity using squeezed input 
light as the optically mediating light field and conclude 
with some general remarks and outlook. 



II. MODEL 

We consider a multimode optomechanical system de- 
scribed in terms of a set of AT modes of effective masses 
rrij , bare frequencies ujj and damping rates jj . The inter- 
action between these modes is assumed to be mediated 
by a single mode of the optical resonator, of nominal 
frequency u c and driven by a monochromatic external 
field of frequency ujl ■ A realization of this system could 
be a high-finesse Fabry-Perot resonator with a series of 
internal non-absorptive dielectric membranes such as de- 
picted in Fig. [TJ a collection of ultracold atomic samples 
trapped at or near the minima of an optical lattice inside 
such a resonator, one or several Bose-Einstein conden- 
sates trapped inside a Fabry-Perot cavity, perhaps with 
a moving end-mirror, or a number of similar setups op- 
erating either in the optical or in the microwave regime. 

The Hamiltonian governing the system is 



,(2) 



linear and 



H = H. 



opt 



H] 



loss ; 



(1) 



where 



H opt = hw c tfA + ih(7 1 e- iuJLt A <! - 7?V" Lt i), (2) 

describes the cavity field mode, pumped externally by a 
field of frequency ujl at a rate rj, 



(3) 



3=1 



describes the mechanical modes, with [Bj, = Sjk- 
The optomechanical interaction is 



H om = -HA^A J2 .90,, (Bj + B)) (4) 



in the linear case and 

H om = htiAY,af}{B j +B]y 



E 



(5) 



in the quadratic case, with g j and g 
quadratic single-photon optomechanical coupling coeffi- 
cients, respectively. Finally, H\ oss describes the inter- 
action of the cavity field and the mechanics with their 
respective reservoirs and accounts for dissipation. For 
simplicity, we will assume that all mechanical modes ei- 
ther couple purely linearly or quadratically. 



III. LINEAR COUPLING 

This section discusses the case of linear coupling in- 
cluding modifications to cold damping and optical spring 
effects. 

With 

A = Ae~ iULt (6) 

and for the optomechanical coupling (0| the Heisenberg- 
Langevin equations of motion in the standard input- 
output formalism [3l| are 



A 



N 

-E 
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A + rj + iA V g , j {Bj +B]) + ^, 



KOi 



Bj = -iuijBj - ^{Bj - B)) + ig ,jAU + ilj, (7) 
where 

A c = uj l - uj c , (8) 

k and 7j are the decay rates of the cavity field and me- 
chanical modes, respectively, and a; n and £j the asso- 
ciated Markovian quantum noise operators. Their two- 
time correlation functions will be given explicitly later 
on. 



A. Linearization 

If the system is driven by a classical field, it is use- 
ful to decompose the various field operators as the sum 
of their expectation values and small quantum fluctua- 
tions, whose dynamics are treated to lowest order only. 
Specifically we expand the operators A and Bj as 



A = (A) + 5 A = a + a, 
Bj = {Bj)+SBj = p j + b j , 



(9a) 
(9b) 



where [a, a^] = 1, [bj, b k ] = djk and we assume that their 
expectation values a and (3j are much larger than the 
amplitudes of the fluctuations, for example, \a\ 2 ^ (a^a). 
(Note that rigorously speaking this assumes all quantum 
contributions to have zero expectation value. However, 
the treatment remains valid for fluctuations with non- 
zero, but small expectation values [Hj].) This yields the 
mean field equations of motion 



a = iAa a 

2 



P j= -iUjPj-?±(Pj-p*)+ig Qt j\c 



(10) 
(11) 
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with steady-state values 



'/ 



where 



-iA + k/2 
2 9o,j, ,2 



A c + $>oj(& + /3*). 



(12) 
(13) 

(14) 



Here we have chosen the phase of r\ such that a s is real 
without loss of generality for the case of cw pumping. 

As is well known, dynamical back-action gives rise to 
frequency shifts - the optical spring effect - and radia- 
tive dam pin g or amplification for each mechanical oscil- 
lator |33|-|35| . In order to capture these effects explicitly, 
we switch to an interaction picture with the transforma- 
tions 



-ifj t 
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where 



Uj = Ldj + Q 



(15) 
(16) 

(17) 



and the expressions for the dynamical shifts fij will be 
determined self-consistently later on (36|. The equations 
of motion for the fluctuation operators then read 



iA- 



M 

a + 1^2,9,^3 

3 = 1 



(18) 



3 2 



bj + ^Ste 2 ^'* + i^e^'V + a) + if 
2 J 



31 

(19) 

where we have introduced the amplified optomechanical 
coupling strengths, 



9j = 9o,j 



B. Elimination of the cavity field 



(20) 



In the physically relevant regime where the optical loss 
is much larger than the mechanical decay rates, k 3> 7j, 
the cavity field follows the dynamics of the mechanical 
oscillators adiabatically. For time scales long compared 
to ft -1 , the formal solution of Eq. (fT8")) reads 



a{t) ~-J2d3 

3=1 



t b- 

u 3 



e lv ^b] 



(A + Uj) + iK/2 (A-Vj) + iK/2 



(21) 



where the noise operator fi n (t) is given by 



/in(*) =yfc dr e^^-^a-^r). (22) 
Jo 



Substituting the formal solution into the equations of mo- 
tion for the mechanical oscillators, Eq. (TT9"]) . gives 



S, + libU 2i ^ 



N 

^9jk 

k=l 
Af 

k=l 



till. — 



where 



9jk = g 3 /gk- 



(23) 



(24) 



and the radiation-induced frequency shifts Qk and damp- 
ing coefficients Tfc are given by 



^fc = 9k 



A-i/fe 



+ 



A + z/ fe 



(A-^.) 2 + K 2 /4 (A + ^) 2 + K 2 /4 



(25a) 



I\ 

2 



k/2 



k/2 



(A + ^ fe ) 2 + K 2 /4 (A — i/ fe ) 2 + k 2 /A 



(25b) 



Due to the appearance of the shifted frequencies z^fe(A) 



in the denominators on its right hand side, Eq. (|25al) is 
a fifth-order polynomial equation for fife (A), so that it is 
not Lorentzian in general. In the strong coupling regime 
gk > k/2, it becomes multivalued for a range of effective 
detunings, a feature that can be interpreted in terms of 
normal mode splitting [f| . In this situation the mechani- 
cal and optical modes are strongly hybridized and it is not 
meaningful to treat them separately. Also, the counter- 
rotating terms in Eqs. (|23[) . combined in particular with a 
driving field blue detuned from the resonator resonance, 
can result in instabilities in the strong coupling regime 
and a concomitant would lead to the break-down of the 
linearization process. 

The strong coupling and blue detuned regimes will be 
the object of future work. Here we focus on the simpler 
case of weak coupling, gk <C k. In this limit the frequency 
shifts fifc remain small compared to the bare frequencies 
of the oscillators, so that we can approximate by LOk 
on the right hand side of Eqs. (|2"5]). see Eq. ([T7|l. The 
frequency shifts fifc are then the sum of two dispersion 
curves while the decay rates Tk consist of two Lorentzians 
with opposite signs, all of which are centered around A = 
±Wfc. The weak coupling regime also allows us to drop 
the counter-rotating terms in Eqs. (|23|) . 

The resulting equations of motion are 



iSlj - ^ 
3 2 



A' 



9jk 



k=l 



iYl k + — 



'</,<"''( A: • fL) ' ^ 



(26) 
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These equations describe a system of N pairwise, linearly 
coupled oscillators, with coherent contributions propor- 
tional to the optically induced frequency shifts fifc, cold 
damping coefficients I\. familiar from the single-mode sit- 
uation, and associated noise operators (/j n + f? n ) and 
£j . The fact that each mechanical element experiences a 
different frequency shift and radiative decay, depending 
on its optomechanical coupling strength and the pump- 
cavity detuning, is the key ingredient that allows one 
to couple oscillators with different bare frequencies on- 
demand, by choosing appropriate pump-powers and de- 
tunings. We discuss selected consequences of this coupled 
dynamics in section V for the case of a two-mode system. 
But first, the next section extends our discussion to the 
case of quadratic optomechanical coupling. 



Zero displacement {fij + /3*) s = is always a solution to 

these equations, and for g^ '■ > 0, which is the case for 
membranes located at minima of the intracavity field, it 
is the only stable solution. However that solution may 

(2) 

become unstable for negative couplings, <?q • < 0, occur- 
ring when oscillators are located at cavity-field maxima. 
In this case, each oscillator settles in one of the stable 
positions, 

±5 (32) 

located symmetrically around the intracavity field max- 
imum. As we will sec in the following sections, these 
two geometries result in qualitative consequences for the 
behavior of the system. 



IV. QUADRATIC COUPLING 



A. Local maxima 



We now consider the quadratic optomechanical cou- 
pling 



.v 



(27) 



3=1 



see Eq. (jSJ), where g$ ■ are single-photon coupling coef- 
ficients. In that case the Heisenberg-Langevin equations 
of motion are, in the frame of the laser frequency u>l, 



A 



Bj = -iUjBj - | (4- - B}) - 2ig$AU(B j + B)) + i| 

with A = Aexp(-iu)Lt) as before. With the expan- 
sion © we have 



- A (2) 

iA y 'a — —a + 

7j 



J 3 — —J i- J 2 
where the effective detuning is now 



(3, ~ -kjjft - f % - - 2i ff g>M a G8 j + /?*), (28) 



3*)\ 



(29) 



and consistently with the linearization procedure we have 
neglected in the second line the vacuum expectation value 
contribution (bw) = 1. The steady-state expectation 
value of the cavity field amplitude is 



'/ 



-iA( 2 ) +k/2' 



(30) 



and the steady-state displacements of the mechanical os- 
cillators are found from the AT coupled fifth-order equa- 
tions 



(4 5 ( 2 j|a s | 2 +^)(/3,+/3;) s = 0. 



(31) 



To simplify the discussion we assume that all coupling 

(2) 

constants g$ '■ have the same sign, and consider first the 
case where they are negative. The stable steady state is 
degenerate, as each oscillator can have a negative or pos- 
itive mean displacement. A particular choice of steady 
state will break the symmetry of the sign of the displace- 
ments and influence the physics of the system. Experi- 
mentally, it could be controlled by addressing the indi- 
vidual oscillators through an additional field or preparing 
them slightly displaced towards the desired steady-state 
position. The Heisenberg-Langevin equations of motion 
for the fluctuations can be approximated to lowest non- 
vanishing order, as 



zA( 2 > - - 
2 



3=1 



111 ■ 



Ka 
iw 



2ig^n r 



2 



Aigyxj(a + at) + 



(33) 



(34) 



where we have introduced the amplified quadratic op- 
tomechanical strengths 



J2) 



(2) 

1 

1,3 



(35) 



and the shifted detuning of the cavity field and mechan- 
ical frequencies are 



A ( 2 )=A( 2 )-4^.^ 

3=1 



(2)^2 
3 3 



A, 



3=1 



jXj, (36) 



(2) - 

Jj + 2g 0j n c 



(37) 



Note that in contrast with the situation for linear op- 
tomechanical coupling the correction Sj^'jLi ffaj^ 2 to 
the detuning A c is now twice as large as the correction to 
the mean- filed dynamics, compare Eq. (|3T|) to Eq. (|29|) . 
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This second-harmonic generation process, a direct acous- 
tic analog of optical second harmonic generation, is a di- 
rect consequence of the nonlinear nature of the quadratic 
optomechanical coupling and is driven by the mean-field 
mechanical oscillations. 

Assuming that the resulting frequencies of the mechan- 
ical oscillators are 



These equations have the same form as in the case of 
linear optomechanical coupling of the mechanical oscil- 
lators, the main new features in this regime are second- 
harmonic generation and the flexibility to choose the sign 

(2) 

and strength of the coupling coefficients g- k through the 
equilibrium positions Xj and x k , see Eq. (|40l) . 



and following the same approach as before the equations 
of motion for the mechanical modes become 



h = [<n 3 - - f 



E4 



*;=i 

JV 



E~ (2) 

k=l 



,9)k 



o rfc 
i\l k — 



(38) 



where the noise operator associated with the adiabati- 
cally eliminated optical field is 



(39) 



and 



B. Local minima 



The situation is changed qualitatively in the case where 
the mechanical oscillators are located at minima of the 
cavity field. The single-photon, single mode quadratic 
optomechanical coupling coefficients are then positive 
and the cavity field results in a tighter confinement 
through the static optical spring effect. The only stable 
steady-state displacement is given by Xj = so that the 
effects described in the previous section are highly sup- 
pressed. On the other hand, a number of novel features 
appear as a result of the fact that the first non- vanishing 
couplings are of higher order in the quantum fluctuations. 
In particular, instead of the second-harmonic generation 
discussed in the previous section, the dominant nonlinear 
effects are akin to four-wave mixing in optics. 



(2) 



(2)- 

(2)- 
9k x k 



(40) 



To lowest order in the fluctuations, the Heisenberg- 
Langevin equations of motion become 



We then find for the radiation-induced frequency shifts 
£l k and damping coefficients Tk 



9., 



- v k 

(A(2) _ Uk y + K 2 /A (A( 2 ) + V k f + K 2 j\ 



A( 2 ) 



(41) 



k/2 



k/2 



(A( 2 ) +lS k ) 2 + k 2 /4 (A( 2 ) -v k ) 2 + n 2 /A 



(42) 



iA^-^a-i^gP^ + bf + 2%) 

3 = 1 



bj ti- 



2igf\a + a v ){b ] +b))+ii J1 



where the effective detuning A^ 2 ) is now 



(44) 
(45) 



(2) — 

In the weak coupling regime, g) Xj -C k, we can ne- 
glect the fast oscillating terms on the right-hand side of 
Eq. (f38|) so that the equations of motion for the mechan- 
ics can be further simplified to 



3=1 



(46) 



3 2 



A/" 



(2) 



fe=l 

j.Vjt(f_ | i\ 



9jk 

4i 5 f ) %e^*(/ in + 4) + ^ 



iO fc + -4- 



(43) 



Applying the same arguments and approximations as be- 
fore to adiabatically eliminate the optical field dynamics, 
we find the coupled equations of motion for the mechan- 
ical modes 
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p e 2iv 3 t _ (2) 

k=l 



l\l k - — 



- V o (2) 

" fc=l 



(e^bfb, + e^^Hfb]) - i £ ^Mblhb 3 + b\b k b)e^) 



k=l 



(47) 



where the mode-coupling is now given by 

(2) (2) / (2) /.on 

9) k =9j J /9 k , ( 48 ) 

the frequency shifts and radiative damping coefficients 
are 



V k = (2g { k 2) ) 2 



(49) 



A( 2 ) - 2v k 



A( 2 ) + 2v k 



(A(2) _ 2v k y + k 2 /4 (A( 2 ) + 2v k ) 2 + k 2 /A 



(50) 



k/2 



k/2 



(A(2) + 2^)2 + K 2/ 4 ( A (2) _ 2^)2 + K 2/ 4 



A fc = (2g 



(2)n2 



2AW 



fe > (A(2))2 + K 2/ 4 ' 



(51) 



/*(*) = dTe^ l2) ^ 2 ^a in (r). (52) 



and 



Note that in Eq. (|47p the quantum noise stemming from 
the elimination of the optical field is now multiplicative. 

Since we can neglect fast rotating terms in the weak 
coupling regime, Eqs. (|47[) further simplify to 



b j = in,- - ^ 



1 M 

iV* ( 2) 

2 5 o,jfc 



A/" 



~ 2 £ & A kb k hbj 

k=l 



iQi- H — — 



- 2i 5 f (/ in + + i>h 2W]t ) + (53) 

These coupled nonlinear equations are indicative of four- 
wave mixing processes, and are formally similar to situ- 
ations driven by x nonlinear susceptibilities in nonlin- 
ear classical and quantum optics (37I . |38| . In that con- 
text, four- wave mixing phenomena are known to lead to 
a rich spectrum of effects that have found numerous ap- 
plications in applied and fundamental studies. They in- 
clude for example sum frequency generation, high har- 
monic generation, optical parametric amplification, self- 
focusing and defocusing, the generation and propaga- 
tion of optical solitons, phase conjugation, the generation 
of quantum mechanical squeezing and entangled photon 



pairs, and more [37|, |38|. The analysis of this section 
paves the way to investigate a similar range of studies in 
quantum optomechanics and quantum acoustics, and as 
such, opens the way to an intriguing and rich direction of 
investigation that will be carried out in follow-up stud- 
ies. One particularly attractive feature of optomechan- 
ics in this context is the ease with which the sign and 
strength of the nonlinear interactions can be adjusted in 
combination with the functionalization of these devices. 
The following sections discuss a few illustrative examples, 
concentrating for now on the simple situation of a two- 
mode system with linear optomechanical coupling only, 
and demonstrating a example of functionalization that 
involves the optically mediated coupling of a nanome- 
chanical system and a quantum degenerate atomic sys- 
tem. 



V. EXAMPLE - TWO-MODE SYSTEM 

Elementary multimode effects can be studied in the 
case where two mechanical modes interact linearly with 
a common cavity field, i.e. Af = 2. In the weak coupling 
regime the Heiscnberg-Langcvin Eqs. (|26|) become 



h = --r e ,i5i - 



«n C; i + -r Ci i 



b 2 e l{vi - v - 2)t 



i92e il/lt (f in + fL)+iii 



1, 



-T P . obi 



in Ci2 



l. 



b\e 



-i(v 1 -v 2 )t 



2 

^- t (/ i „ + / i r n )+*6 (54) 
where we have introduced the effective decay rate 



r — r 



7j, 



the cross-damping rate 

r c ,j = 9jk^k, 
and the cross-coupling frequency 

^c,j = 9jk^k- 



(55) 



(56) 



(57) 



The cross-damping rate is associated with effects such 
as sympathetic cooling or heating of coupled oscillators, 
while the cross-coupling frequency is associated with the 
coherent aspects of multimode coupling, such as e.g. 
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quantum state transfer. Note that when the two effec- 
tive mechanical frequencies are matched, v\ = ^2, the 
cross-coupling frequencies and damping coefficients be- 
come equal Cj i = ri Cj 2 , r c .i = r Cj 2. 

In terms of the motional quadratures 



vi, v 2 (a.u.) 



1 

1 



Pi = 



Eqs. ([55)1 become 



u = Mu + S, 



(58) 
(59) 

(60) 



where u = (X\ 1 Pi, X2, A) i the drift matrix M. is given 
by 



-r e ,i 


2i/i 


-r c ,i 


2f2 Ci i 




-r e ,i 


— 2fi Ci i 


-r c ,i 


-r c>2 




-r e ,2 


2^2 




-r c , 2 


-2^2 


-r e , 2 



M = - 



and the noise operator matrix S is 

t ° . \ 


V6e- W2t +.92-Fi„/ 
with the effective optical noise operator 

F in = V« [ dr e^-^^-^Mr) + adj. 
Jo 



(61) 



(62) 



(63) 



In the following we limit ourselves to the case of Gaussian 
states and assume that the mechanical elements are cou- 
pled to statistically independent Markovian heat baths 
at low temperature T characterized by the two-time cor- 
relations functions [33 



7 3 -(2n thJ + l)<5(i-t'), 



(64) 



the 



where fith,j = [exp(huij /fc^T) — 1] . From Eqs. 
first moments of all position and momentum quadratures 
are always zero, and one readily derives a closed set of 
differential equations for the second moments of the os- 
cillator quadratures. The equation of motion for the co- 
variance matrix Vij = ^(uiUj + UjUj) can be written as 



V = MV + VM +D, 



(65) 



where Ai is the drift matrix, and the inhomogeneous 
term D is given by Dj 
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This system of equations can be solved exactly in prin- 
ciple but the resulting solutions lack transparency. More 
physical insight into the coupled dynamics of the oscilla- 
tor can be gained in the resolved sideband limit and in 
the Doppler limit, two cases to which we now turn. 




•40 - 20 
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FIG. 2. Effective frequencies of the first (red, solid) and sec- 
ond (blue, dashed) mechanical modes in the resolved side- 
band regime wi, u?2 S> ft for g\ > g-z- The dotted lines 
denote the bare frequencies of the mechanical modes. The 
shifted frequencies are matched at the intersections. The in- 
set gives a detailed view on the intersections occurring in the 
red-detuned side. The parameters used are wi = 2n x 20MHz, 
lu 2 = 2tt x 19.95MHz, k = 2tt x 1MHz, g 1 = 2tt x 0.3MHz, 
02 = 2tt x 0.12MHz. 



A. Resolved sideband regime 

The resolved sideband regime is characterized by slow 
optical resonator damping compared to the mechanical 
frequencies uj\, L02 S> k. (We assume without loss of 
generality that uj\ > L02 in the following.) It is known 
from the standard theory of single-mode cavity optomc- 
chanics that in the absence of coupling each oscillator ex- 
periences a frequency shift f^(A) comprised of the sum 
of two dispersion curves centered around A = ±Wfc with 
characteristic widths k, see Fig. [21 while their damping 
coefficients r^(A) are the sum of a positive and negative 
Lorentzian, also of width k, and centered at these same 
detunings [ill ]. 

The red and blue detunings for which the two oscilla- 
tors are brought into resonance with each other are given 
by the intersections of the two curves in Fig. [5] There are 
two such intersections on the red-detuned side of the op- 
tical resonance, one near A = —ujk and one on the wing 
of the dispersion curve. Two more intersections are lo- 
cated similarly on the blue-detuned side. At these points 
we have v\ = V2 = v and therefore Ci i = f2 Ci 2 = fi c and 
r c ,i = r c ,2 = r c from Eq. ([57]) . It is for those detunings 
that the coupling between the two mechanical oscillators 
is most effective. 

The detailed dynamics of resonant mode coupling de- 
pends on whether the resonance condition v\ = V2 occurs 
near the center or on the wing of the dispersion curves 
that comprise fJfe(A). Near the resonance A = ±Wfc opti- 
cal damping (or amplification) dominates over the optical 
spring effect, see Eqs. (|25a[) and (|25b[) . and to first ap- 
proximation one can safely neglect that effect compared 
to those of r Cj fc by setting fl c ,k = 0, see Eqs. ([5"6"[) and 
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(|57p. The drift matrix Ai reduces to 



M 



-r e ,i 


2i/ 


-r c 





-2v 


-r e ,i 





-r 


-r c 





-r e , 2 


2;/ 





-r c 


-2v 


-r e 



(66) 



and the dominant consequence of the effective optome- 
chanical coupling between the oscillators is the modifica- 
tion of their rate of cold damping or amplification com- 
pared to their uncoupled values. One such example is 
shown in Fig. [31 which plots the evolution of the variances 
of the motional quadratures of the mechanical modes 
when their effective frequencies are matched around the 
center of the Lorentzian curve on the red side of the cav- 
ity field. The mechanical modes are both assumed to be 
initially in thermal equilibrium. In the absence of cross 
interaction, both oscillators are subject to cold damp- 
ing (dotted and dot-dashed lines) with cooling rates that 
differ since g\ ^ g 2 - The coupling changes the situa- 
tion significantly, with both cooling rates now markedly 
slower (solid and dashed lines) . This is a result of the ex- 
change of thermal excitations between the two oscillators 
that inhibits their direct optical cooling. (Note that our 
numerical calculations suggest that the long-time limit 
and thus the cooling limit does not change.) Performing 
a normal-mode analysis on the two modes reveals that 
the two out-of-phase oscillations decouple from the cav- 
ity field, while the two in-phase modes are cooled with 
an increased damping rate. 

We now turn to the resonances v\ = v 2 that occur 
on the wings of the dispersion curves. In these cases 
the optical spring effect dominates over cold damping 
(or amplification), a consequence of the faster decay of 
the Lorentzian Tfc(A) (with 1/A 2 ) compared to the dis- 
persive nature of Qfc(A), which decreases as 1/A. That 
situation is therefore characterized by the dominance of 
coherent exchange between the two oscillators, and we 
approximate the corresponding drift matrix by 



M = 





V 







— v 


-71/2 


-f2 c 








n r 


-72/2 


V 


-n c 





—v 


-12/2 



(67) 



In this regime, and provided that the noise sources as- 
sociated with the optical field can be managed, one can 
realize coherent effects such as quantum state transfer 
between the two mechanical modes. We return to this 
effect in the Doppler regime, which we consider in the 
following. 



B. Doppler regime 

In the Doppler regime, uj\, u>2 ^ k, the width of the 
optical cavity is so large that it washes out the two dis- 
tinct features of fife (A) and Tfc(A), characteristic of the 
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FIG. 3. Upper plot: short time (linear time scale) and lower 
plot: long time (log time scale) evolution of the motional 
quadratures of the mechanical modes in the single and two- 
mode scenarios. The mechanical modes of frequencies cji (red, 
solid) and U12 (blue, dashed) are both cooled down while inter- 
acting with each other. For comparison the uncoupled cold 
damping of the modes wi (orange, dotted) and LJ2 (green, 
dotdashed) are also shown. Two-mode coupling results in 
slowing down of the individual cold damping rates. Here 
wi = 2tt x 20MHz, uj 2 = 2iv x 19.99MHz, k = 2tt x 0.95MHz, 
gi = 2tt x 50kHz, g 2 = 2ir x 10kHz. 



resolved sideband regime, and one can approximate 



2A 



A 2 + K 2 /4' 

2 2 Akja, 

1j (A 2 + K 2 /4) 2 ' 



(68) 
(69) 



The qualitative behavior of the effective frequencies 
j (j = 1, 2) is plotted in Fig.@]for gi > 32 and shows the 



v 



intersections of the two curves on the red-detuned side 
of the cavity resonance. (They would lie on the blue- 
detuned side for g 2 > gi). Since decays with the 
third power of the cavity linewidth while f2j(A) scales as 
1/k 2 , it follows that for the large cavity damping rates, 
characteristic of the Doppler regime, we can neglect in 
first approximation the effects of cold damping (or am- 
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FIG. 4. Qualitative behavior of the effective mechanical fre- 
quencies for the first (red, solid) and second (blue, dashed) 
oscillators in the Doppler regime. The shifted frequencies are 
matched at the intersections. Here we chose gi > gi. A 
set of parameters is wi = 2n x 100kHz, u)2 = 2n x 93kHz, 
K = 2tt x 1MHz, gi = 2tt x 90kHz, 52 = 2tt x 50kHz. 



plification) and at resonance we find again 
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-71/2 


7' 







—v 


-7l/2 
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n c 


-72/2 


;/ 


-fie 







-72/2 



(70) 



Note that the argument leading to this result is different 
in the resolved sideband limit, where the radiation in- 
duced damping can be neglected if the effective detuning 
is outside of the mechanical sideband, and in the Doppler 
regime, in which case it is negligible due to the large cav- 
ity lincwidth. 



C. State transfer 

The previous discussion showed that both, the resolved 
sideband and the Doppler regimes, provide situations 
where cold damping can be neglected and the two os- 
cillators are resonantly coupled by the effective coherent 
interaction 



H e = nn c (bibl + b[b 2 ) 



(71) 



This is the familiar beam splitter Hamiltonian, which 
is known in particular to lead to quantum state swap- 
ping between the two oscillators. However, an important 
issue is that the optical field mediating the interaction 
also results in the appearance of quantum noise, further 
limiting the fidelity of state transfer in addition to the 
clamping noise associated with the mechanical coupling 
to thermal reservoirs. Clamping noise can be reduced in 
principle in cryogenic environments or perhaps in levi- 
tated structures |42| . What we show in this brief section 
is that the fundamental quantum noise associated with 



the optical field can also be reduced in principle by the 
use of squeezed input fields with correlation functions [43[ 

(a in (t)a in (t')> = e-^ t+t 'hl5(t - t'), (72) 

(^(t)al(t'))=e i ^+^M*S(t-t'), (73) 

(^(t)a Ul (t'))=N6(t-t'), (74) 

(a in (t)al(t')) = (N+l)5(t-t'), (75) 

where ujl is the central frequency of a squeezing device, 
and positive-valued N and complex-valued M are the 
squeezing parameters which define the strength of the 
squeezing as well as its direction in phase space, given by 
the complex phase of M. In the following, we assume an 
ideal squeezed state which satisfies 



M = ^N(N + Tje 



-2i8 



(76) 



The equations of motion for vacuum input noise are re- 
covered by substituting N — > 0, M — s- 0. 

For concreteness we consider the explicit situation of 
a hybrid two-mode mechanical system where one of the 
modes is a recoil-induced side-mode of a Bose-Einstein 
condensate [|[ trapped inside an optical resonator with 
an oscillating end-mirror, and the other is the center-of- 
mass mode of oscillation of that mirror. As discussed in 
Ref. 
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in the case where a system is prepared in such 
a way that the shifted optomechanical frequencies of the 
two oscillators v\ and v 2 are equal, f2 Cj i = f2 Ci 2 = Q c , the 
equations of motion for the quadratures of the coupled 
mechanical oscillators reduce to 
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where the subscripts 1 and 2 refer to the oscillating mirror 
and the condensate side-mode respectively, and the cross- 
coupling frequency is 



2A 



3i52- 



k 2 /4' 



(77) 



For a squeezed input field the two-time correlation func- 
tion of the effective noise can be approximated as 

(An(t)An(*')+An(t')An(*)) K 



A 2 + k 2 /4 



|M| cos[2(6» ;j - e c )] +N+-)5(t- t'), (78) 



with 9 C = arg(iA + k/2). Here we have assumed that 
intrinsic mechanical decoherence is negligible for the time 
scales of interest. 

Figure \5\ shows the co-evolution of the variance of one 
of the motional quadratures of the two oscillators. In this 
example the resonator is driven by a coherent field and 
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FIG. 5. Motional quadrature of the moving mirror (red, solid) 
and the BEC (blue, dashed) in the presence of vacuum noise 
for wi = 2tt x 101kHz, u 2 = 2?r x 100kHz, k = 2% x 1MHz, 
t/i = 2tt x 100kHz, p 2 = 2tt x 10kHz. 

<AX 2 > 
1.0 r 




FIG. 6. Motional quadrature of the moving mirror (red, solid) 
and the BEC (blue, dashed) for squeezed vacuum noise with 
N — 1 and squeezing phase 8 S — 9 C — n/2. Same parameters 
as in Fig [5] 

the initial states of the BEC and the moving mirror are 
the ground state and a thermal state with mean phonon 
number n — 1, respectively. As expected, the fidelity 
of the coherent quantum state transfer, which would be 
unity at times (it/2 + pTr)Cl~ , with p integer, in the ab- 
sence of noise, is increasingly reduced by the quantum 
fluctuations of the light field, even without intrinsic me- 
chanical decoherence mechanisms. 

The use of a squeezed optical field to reduce the detri- 
mental effect of quantum noise is illustrated in Fig. M 
where the improved fidelity is readily apparent in partic- 
ular at time t — 7r/2fi c . As illustrated in Figs. [7] and [51 
the state transfer fidelity depends strongly on the phase 
difference between the cavity field and a squeezing device. 
For 6 S = c ±n/2, the fluctuations in the position quadra- 
ture of the cavity field are minimized so that very high 
fidelity state transfer is achieved, T = 0.997 for N = 10 
in our example. As would be intuitively expected, we also 
find that strong squeezing leads to an increased sensitiv- 




FIG. 7. State transfer fidelity for squeezed vacuum noise with 
N — 1 (red, solid), N = 10 (blue, dashed) and vacuum noise 
(gray, dotted) as a function of the phase difference between 
the squeezing input and the cavity field. The initial states of 
the BEC and the moving mirror are the ground state and a 
thermal state with mean phonon number n — 1, respectively. 



F 




FIG. 8. Same as Fig. [3 but for the BEC initially in a squeezed 
state with position quadrature variance of 0.025 and momen- 
tum quadrature variance of 10 and and the moving mirror in a 
thermal state with mean phonon number n — 1, respectively. 



ity of the state transfer fidelity on that relative phase. 
As a comparison to Fig. which is for a BEC side-mode 
initially in the ground state and the oscillating mirror in 
a thermal state with n = 1, Fig. [8] displays the fidelity 
of state transfer where the BEC is initially in a squeezed 
state with (AAf) = 0.025 and (AP 2 2 ) = 10. The transfer 
fidelity in that case exhibits an increased dependence of 
the relative phase between the cavity field and a squeez- 
ing device, as would be expected. In the optimal case 
6 S = 9 C ± ir/2 and N = 10, the fidelity of squeezed state 
transfer is T = 0.978, which is somewhat lower than for 
a vacuum state, a consequence of the higher sensitivity 
of squeezed states to decoherence. 
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VI. CONCLUSION 

In summary, we have analyzed the optically mediated 
interaction between multimodc phonon fields in cavity 
optomechanical systems where it is possible to adiabat- 
ically eliminate the electromagnetic field. A number 
of possible forms of this interaction were identified, in- 
cluding multimodc beam-splitter interactions and, in the 
case of quadratic optomechanical coupling, a plethora of 
possible nonlinear effective interactions, such as analogs 
of second-harmonic generation and four-wave mixing in 
nonlinear optics. The specific forms of the effective in- 
teraction depend sensitively on the sign of the couplings. 
Importantly, we found that incoherent effects due to 
quantum fluctuations - associated primarily with the 
elimination of the electromagnetic field - need not over- 
whelm coherent mode-coupling effects. As a result one 
can expect that analogs of nonlinear effects such as four- 
wave mixing, soliton generation, the creation of entan- 
gled phonon pairs, and others should be achievable in 
quantum acoustics. 

Our treatment was kept sufficiently generic that it can 
be applied to a wide variety of systems and parameter 
regimes, and as a concrete example we have discussed 



quantum state transfer between a momentum side-mode 
of a Bosc-Einstein condensate and the center-of-mass 
mode of oscillation of the end mirror of a Fabry-Perot in- 
terferometer. We showed how decoherence effects can be 
reduced significantly in that case by the use of squeezed 
optical fields. 

Future work will expand these results to more general 
situations and apply it to the study of practical situations 
where the effective mode coupling is nonlinear and inves- 
tigate particular effects such as phase conjugation, en- 
tanglement, and the creation of entangled phonon pairs, 
for example. We will also exploit pulsed optomechanics 
ideas to further expand the toolbox of quantum optome- 
chanics. 
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